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1. Introduction 

A/" = 4 supersymmetric Yang-Mills theory attracts a lot of interest after the discovery 
of AdS/CFT correspondence. During the early years, due to the BPS properties of the 
states from dimensional reduction of the supergravity, studies on the CFT side were mainly 
focused on the protected operators (see [1] for a review). The proposal of BMN correspon- 
dence [2] provides further evidence of AdS/CFT beyond supergravity. Operators Tr(Z'^) 
are BPS and correspond to the string vacuum states |0,p^) in the light-cone quantization 
of the string in a pp wave background, whereas operators like Tr((/)jZP0jZ'^~P) are not 
protected and get mixed in order to remain the eigenstates of the one-loop dilatation op- 
erator. In the BMN limit these eigenstate operators (BMN operators) correspond to the 
string excitation states anOpn\Q^p'^) (see e.g. [3] for a pedagogical review). BMN operators 
involving vector [4-6] and fermions [7-10] are also examined. Exact BMN operators at 
finite N and J are obtained in [11] and extended to the whole superconformal multiplet 
incorporating BMN operators with various combinations of two impurities by using su- 
perconformal symmetry. Operator mixings beyond one loop in the scalar sector are also 
discussed in [12]. 

Unlike the mixings in the pure scalar sector, mixings between scalar and fermionic sec- 
tors are less studied. In [7], some low dimension operator mixings are discussed. Instanton 
contributions to their mixings are also studied in [13]. After mapping the operators into 
spin chains and the dilatation operator into a spin-chain Hamiltonion, a Bethe ansatz can 
be used to solve the mixing problem [14,15]. Integrability and the analytic Bethe ansatz 
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make it possible to obtain anomalous dimensions for operators with more than two impu- 
rities. For = 4 SUSY YM, this was first proposed in [14] for an so{6) spin chain and 
further extended to the psu{2,2\4) super spin chain in [16]. A subsector su(2|3) dynamic 
spin chain which includes the mixing between fermion pairs and bosons is discussed along 
these lines in [17]. The coordinate-space Bethe ansatz in calculating the spin-chain S- 
matrix can also produce the eigenstates with multiple impurities [15,18]. This is known as 
the perturbative asymptotic Bethe ansatz in this context. The sm(2|2) dynamic S-matrix 
of the su{2\3) subsector was subsequently discussed using the nested Bethe ansatz [19]. In 
this subsector, the fermions form singlets of the SU{3) subgroup of the R symmetry. 

In the present paper, the fermions we will discuss transform under the 5*0 (4) = 
SU{2) X SU{2) subgroup of the SU{4)Fi symmetry as two fundamental representations 
for the two 5'[/(2)'s respectively. As a first attempt, we only discuss the mixings between 
BMN operators with two impurities at planar level and the fermion pair form a Lorentz 
scalar. In particular, we will discuss the mixings between operators: 

Or'''= (^)'^'Tr(ZZ^+i) (1.1) 



O^^f = (-^^ TT{(t>.ZP<P,p-n , for i,j = 1,2,4,5 (1.2) 

^Pf=(^^y(^^^ T¥(V'^'"ZPV^„Z^-P) , fori,j = l,2, (1.3) 

= (^^^ Trii;,,^ZPi;fz'-n , for z,j =3,4. (1.4) 

We have normalized the operators so that 



for tree level planar diagrams. From the SO (4) point of view, since scalars with two 
impurities can form single trace operators of 1 + 3 + 3 + 9 representations, operators like 
( |1.3| ) with i = 1, 2 and j = 3, 4 which is in vector 4 representation can not be mixed with 
( |1.2| ). Because the common representations of two scalars or two spinors are 1 + 3 + 3, we 
expect that they only mix in these representations. 

In [20], a method is developed to obtain the mixing between these operators using the 
property that the eigenstate operator should be annihilated by some supercharges. We 
will reproduce their results by directly diagonalizing the dilatation operator and give all 
the mixing patterns. Since we only consider two impurity cases, it is easy to do it without 
using Bethe ansatz. We will follow the method used in [21-25] to deduce the dilatation 
operator in the name of effective vertices. Even though the algebraic method was used to 
extend the scalar one-loop dilatation operator to the full AA = 4 SUSY YM [27] , we will still 
deduce the one-loop effective vertices by feynman diagrams since it is not hard to calculate 
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the correlation functions with only scalar fermion pairs [9, 10] and we can formulate the 
effective vertices taking advantage of the scalar property of the fermion pair. We also need 
to calculate the effective vertices to 0{g^) which is the lowest order at which the mixing 
between scalar sector and fermionic sector occurs. 

In section ^, we review the matrix model effective vertices and their relation with 
dilatation operators and calculate the one-loop effective vertices with one scalar fermion 
pair. In section ^, we review the mixing of pure bosonic operators using the one-loop 
effective vertices and obtain the mixing of operators with one scalar fermion pair using the 
results from pure scalar operator mixing. In section we perform a general order analysis 
about the form of the eigenstates. In section |5|, we give the matrix model vertices to 0{g^) 
with a pair of fermions. In section |6|, we use the effective vertices to obtain the dilatation 
matrix to 0{g^) and solve the eigenvector operators. The results are consistent with [20]. 
Section |^ is the discussion and outlook. Appendix ^ gives the conventions we use in our 
paper. 

2. One-loop matrix model effective vertices and dilatation operators 

Let us recall the relation between matrix model vertices and dilatation operators in [25]. 
The one-loop two-point function can be written as 

(O„(x)O^(0)) = ^^i^(5,^ + T„^log(|xA|-2)) . (2.1) 

The dilatation operator matrix element can be expressed as 

= ( J + 2)5/ + T^^{S-y^ . (2.2) 

The basic idea of [22] is to construct matrix model effective vertices H to reproduce the 
one-loop T matrix elements for scalar sector by 

= (0„ H O^) = H^^O^Of^) = H^^S^i^ (2.3) 

and 

H = ^ {Vn + Vp + Vk) (2.4) 

with 

Vd = 1- </>-][</>+, (!>-] ■■ +N : Tr0+,/.- : - : T:(,/.+ )Tr(</)- ) : (2.5) 

Vf = Tr[0+,(/>+][0-,(/>~] : (2.6) 

Vk = -\:Tv[cP+,a(Pl,^-]: (2.7) 

where repeated m and n sum over 1, . . . , 6. and denote the (p at x and respectively 
and contractions are nonzero only between 0+ and (j)^: 

{{^t)ab{(p])cd) = S^,5ac^M , (0+^) = (0+0+) = . (2.8) 
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It can be proven that Vd is cancehed out when acting on traces of scalars (j). But we need to 
reexamine this result when considering operators with fermion pairs hke later. From 



( |2.1| ) , we see that Ha'^ are the matrix elements of one-loop Dilation operator when acting on 
the operators with traces of scalars. For the operators ^'fj'^ with a scalar fermion pair the 
forgoing discussion does not need to be modified except that the H operator should include 
the fermionic sector. Notice that for states in (1.4), S^^ is diagonahzed and normalized. 
So to change the operators into Matrix model operators we need to make the replcement: 



(2.9) 



Since we are interested in the operators with only scalar fermion pairs, we need only to 
calculate the correlation functions with scalar fermion pairs in order to obtain the effective 
vertices. The matrices in the fermion propagators of SUSY YM would be assembled into 
matrix model coefficients and will not appear in the correlation functions of two fermions 
in the matrix model. We can specify the contraction rules of the fermions in the matrix 
model: 

{{^a)ab{i^a)cd) = eaaSad^bc (2-10) 

and require 

11^1 = r, {ry = -i^a. (2.11) 

The correspondence between the SUSY YM operators and matrix model operators is: 

(2^/^-)^-(^)^. (2.12) 
We need to calculate the one- loop diagrams in SUSY YM in figure |^ to get the effective 



1^ 



Figure 1: One-loop diagrams in calculation of the effective vertices. The scalars and fermions are 
denoted as solid and dashed lines respectively. 

vertices. One can consult [10] for a detailed calculation of these diagrams. Because we are 
considering operators with only one scalar fermion pair, we do not need to worry about 
the contractions between different scalar fermion pairs. For a connected diagram like: 



4^g^ ,:.l.^0a ix-y),^H-ix-y) ^) 



we can add another free fermion propagator to combine with the fermion su{2) indices to 
form a Lorentz scalar coefficient. As a result we include another trace of a fermion pair in 
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the matrix model vertices. The coefhcients are obtained so that the T matrix is the same 
both from SUSY YM side and the matrix model side: 



H. 
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2(47r2) 



-- : Tr ({^^^-, }[Z+, Z-]) TVV^^^+ : + 
+ 1 : TY {[i;^^-,Z+][Z-,^P^^]) TrV^f-V^J^ : + 



■t=l,2 



N : ( Tr'0"-V^+ - -^TrV'""TrV'^+ J Tr'0;;V'J : +(+ ^ 



+ ^ 

+ ^ 

+ ^ 

+ ^ 

+ ^ 

+ ^ 

+ ^ 



-i ^ : TV ([V''^^+,Z-][Z+,^r]) TrV^f+Vi- 

i=3,4 

+ i ^ : TV {[^l;^^+,Z+][Z-,^t]) ^V'^V^^ 



+ 



+ 



i=l,2 



- ^ : TV ([i^P'-,Z-][Z+,{4-']) TrV'f+Vij^ : + 



i=3,4 



+ ^:TV({^/^^+,V;^}{<+,V5t^}) : + 



1 



where (+ — ) means exchanging the + and — in the superscripts. We will use this 
notation throughout the paper. 

Acting on an operator ^'^j , the terms in (^^) which contribute are only in ( |2.5| ) and 
can be written as: 



Vd ~ ^Ti {[Z^ , Z-][Z+ , Z-]) - NTtZ-Z+ + TtZ-TiZ^ 



(2.15) 



Like in the bosonic case, we can prove that these terms cancel with the second and 
third line in ( 2.14j ). To see this, contract one Z~^ in the first term of (2.15) with Z in 
Tr{ij''°'ZPipiZ-^~P) and after some cancellations, the terms left are 

i |-TV(V''"-[Z+, Z-]ZPi;i- Z-^-P) + ^Tr(V''"-Zf 

-Tiii;'''- ZPi;i-[Z+ , Z-]Z'^-P)] + ^Tr{i;'''- ZPtP^- Z-^-P[Z+ , Z'])^ 

= -iTV([7/.*"-, [Z+,Z-]]ZPiji-ZJ-P) - ^Tt{^'''- ZP[^i- , [Z+,Z-]]ZJ~P) . (2.16) 

Contracting the first term in the second line of ( p. 14 ) with the same operator and keeping 
track of the minus sign when one fermionic operator moves over the other, one gets 

^Tr([V^--, Z-]]ZPi;i-Z'-P) + iTr(V^— Z^^", Z-]]Z'-p) , (2.17) 

which just cancels ( |2.16| ). Also as in the bosonic case, contracting the Z^ which is used 
in the contraction in ( ^.161) with the Z~ inside the first term in ( 2.15| ) cancels the other 
terms left. 
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3. One-loop operator mixing 



We only need to discuss the operator mixing at planar level. The pure bosonic operator 
eigenstates are given in [11]. Since wc confine ourselves to only two impurity cases, it is 
easier to directly diagonalize the H matrix than use a Bethe ansatz. Let us review the 
bosonic operator mixing first. Oij can be decomposed into 50(4) representations: 

o: 



= -Tr(ZZ-^+^) , 



1=1,2,4,5 



1 I o^r^ + O^! 



3,J 



-Of" 
2 ^ 



OP:i = - ( QP'-^ - QP'A 



(3.1) 



^'-^ and Of-^ are 50(4) singlets. O^'/^ is in the 9 representation and O^^^^ is in the 3 + 3 



(ij) 

representation. 

The H matrix for Of"^, {p = -1, 0, 1, 



,[J/2]) is 



for J odd: 

/ 1 -1 
-1 2 -1 

g^N -1 2 -1 
47?" 

-1 2 -1 

V -11/ 

Eigenstates: n = 0, • • • , (J + l)/2 



(^2cos j^,..., 2cos 



(2p+3)7rn 
J+3 ' 



...,2cos(^ 



for J even: 
/ 1 -1 
-1 2 -1 
p2jv -1 2 -1 

47r^ 

-1 2 -2 
V -12/ 
Eigenstates: n = 0, • • • , J/2 + 1 
(2cosj^,... ,2cos(^,..., 

r, (J+l)7m \ 

2 cos ^ j_^_3 , cos 7rn 1 



(3.2) 



where all the unspecified matrix elements are zero. The eigenstates can be written in both 
cases as 



0/'" = Vcos (2^ + 3)^" 0^ + 2 



p=0 



cos 



vrn 
J + 3 



The eigenvalue is 



6Do 



9^. ■ 2 
^ -4 sm^ 



(3.3) 



(3.4) 



47r2 J + 3 

In fact, we could take the coefficients before Of '"^ to be cos or sin which 

satisfies the eigenvalue requirement except the last line and the first line in the matrix 

for J = odd. The last line and the first line act as boundary conditions: the last line 

:p+a)n 
J+a 



determines the coefficients to be cos if require a, b to be independent of J, 



the first line determines a = 3. 
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For the symmetric operator Of^ij-^, the H matrix is the same form as in the Oi case 
except that p starts from not from —1. As a result, the eigenstates are: 



Cj,n _^ (2ff + l)7rn J 



(3.5) 



and the eigenvalue is 



p=0 



SD2 = sm^ 



47r2 J + 1 

For the anti-symmetric operators O^.^-^, {p = 0,1, . . . , [(J — l)/2]), the H matrix is 



(3.6) 



for J odd: 

/ 2 -1 
1 2 -1 

a^N -1 2 -1 



V 



-1 2 -1 

-1 3 / 

Eigenstates: n = 1, • • • , (J + l)/2 
(2sin^,..., 2sin(^, 



for J even: 
/ 2 -1 
-1 2 -1 

g^N -1 2 -1 



-1 2 -1 
V -12/ 
Eigenstates: n = 1, • • • , J/2 , 



(3.7) 



(^2 sin 



27rn 
J+2' 



, 2sin 



{J+l)TTn 
J+2 , 



^ . (2p+2)7rn 

...,2sin^)^ 



The eigenstates can be written as 



^sin 
p=0 



{2p + 2)7rn j 



J + 2 



and the eigenvalue is 



(5£>2 



4sin^ 



nvr 



(3.8) 



(3.9) 



47r2 J + 2 

For operators with a scalar fermion pair, it can be checked that the Hamiltonian can 
be recast as 



R 



2(47r2) 



(3.10) 



i=l 



where Hi^j+i is the graded permutation operator which exchanges the adjacent fields and 
picks up a minus sign if the two fields are fermions. L is the total number of fields in the 
operator to be acted on. This equation may not be true for two or more pairs of fermions, 
because there could be changes like ijj'^ipaipl^ipg 'ijj°'ijj^ipgipa- 

As in the pure bosonic cases, we can also define symmetric and antisymmetric combi- 
nations for = 1,2) or = 3,4): 



2^ 

2^ »j 



(3.11) 
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^^Ij^ are in the 3 representation for i,j = 1, 2 or 3 representation for i,j = 3, 4. We expect 
that Oj^l will be mixed with ^'j/)"^ and Of with "^^^ij]'^. 

The H matrix for ^^Ij-j is the same as in Oj'^j^ case, therefore the eigenstates could be 
obtained just by replacing O with ^ in ( |3.5| ) and the eigenvalues are the same. In a similar 
way, the ^'j'^'jj case is the same as Cj^y| case. 

From these diagonalization procedures, we see that for a matrix with —1,2, —1 at the 



near diagonal and diagonal positions like in (|3.2| ) and (3/7), we can use an ansatz like 
cos ^^P+^'l^"''^ or sin ^^^j^^^"^ as coefficients because they automatically satisfy the eigenvalue 



equations. Combining the coefficient before the matrix (|3.2| ) and (3^), the eigenvalue 
is ^^4sin^-j^. The variables a, b and sin or cos are determined from the boundary 
condition: the first line and the last line of the matrix and the requirement that a and 
b are independent of J. This is much like solving an eigenvalue problem of a differential 
equation with boundary conditions. This insight is helpful when we solve the mixing 
between pure scalar operators and operators with fermions. 

4. Order analysis 

In general, the SO{A) singlet part of Of and ^'|'./|~^ should also be mixed with operator 

yP^^ = Tr{D^ZZPD''ZZ-^-P-^) . (4.1) 

In paper [26], the authors use the method developed in [20] to discuss the mixings among 
these operators in which V^'"^ appears only from 0{g^) terms. We will try to understand 
the general form of the eigenstate vector by order analysis in this section. 

Let us look at the Dilatation matrix structure of these three kinds of operator system 
to 0{g'^). For the combined operator basis (. . . , O^'^^, YP'^'"^, . . . ), 
the correction to the dilatation operator up to 0{g^) in matrix blocks is in this form: 



5D 



/0{g^)+0{g^) 0{g^) 0{g^) 

0{g^) 0{g^) + 0{g^) 0{g^) |. (4.2) 
\ 0{g^) 0{g^) 0{g^) + 0{g^~^ 



In general, the eigenvector should look like this 

\OP^^ + Bp^""^-^ + CpVP"^-^ . (4.3) 

According to general perturbation theory of quantum mechanics, we first diagonalize the 
0{g^) part of the matrix. This is done in previous section and in paper [11]. The three 
groups of eigenvectors correspond to / 0, = Cp = or Sp / 0, Ap = Cp = or 
Cp 7^ 0, Ap = Sp = and the eigenvalues are not degenerate. We will use O'^'", 
and y-'~2,n {;ig]2ote these eigenstates. So the non-degenerate perturbation theory could 
be used to calculate the correction to the next order 0{g^). It is easy to see that there is 
no correction to the eigenvalues to this order since the correction terms are proportional 
to (0"^'"|i?(0(53))|0-^:«), (§-^'"|i7(0(53))|§"^:") or (y-^'"|if(0(53))|y"^:«). If we are only 
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interested in the mixing between O^'"^ and ^P'-'"^, we can just take the upper left 2x2 
matrix and diagonahze it. This is what we wih do later. We will see that it is easier to 
do it without using the quantum mechanics perturbation theory result. We can then do 
the perturbation to 0{g'^). Just like the 0{g^) case, the nondiagonal matrix blocks do not 
contribute to the eigenvalue but contribute to the mixing coefficients. This means that if 
we only consider the lowest order of the mixing coefficients between O^'^ and V^''^~'^ we 
can just use the diagonal 0{g^) blocks and nondiagonal 0{g'^) blocks. Combining 0{g^) 
and 0{g^) results, the three kinds of solutions for the coefficients should be Ap ~ 0(1), 
0{g), Cp ~ 0{g^) or ~ 0{g), Bp ~ 0(1), Cp ~ Oig) or Cp ~ 0{g^), Ap ~ 0(5), 
0(1). Notice that the finite J condition plays an important role here, because in 
the infinite J limit the three groups of the leading order eigenstates degenerate and the 
non-degenerate perturbation argument can not be used here. Just diagonalizing the upper 
left 2x2 block matrix would not be valid. 



Bp ~ 
Bp - 



5. Effective vertices of 0{g^) 

In order to study the mixing between Ofj'^ and ^'fj'^ ^ we need to calculate dilatation 
operator changing fermions into bosons and vice versa. The lowest order nonzero diagrams 
are of 0{g^): 



I I 
I I 
I I 



(A) 



(B) 



It is easy to calculate the divergent parts of these two diagrams using relation (B.l) in [20] 
and differential regularization [28]. The effective vertices can be obtained in the same spirit 
as in section ^: 

-.3 



+(+ 

9' 



Tr 



BC 



<— > 

,3 



(5.1) 



rTr 



8(2^)3 
+(+ - -) 

Using ([A. 3D , we can change $ into Z and Z: 

3 r 1 



(5.2) 



9 



" (27r)3 
1 
4 



Tr 



i,[r'-,z^^]{-W+,€-],z^] - [[zt,z'^u'-] + [[v^i-,z+],z'=+]) 



+-\r'-,zi\{[[z^^,i^i-],zt] - [[z^^.ztu^]) - \i^wr'- .ztwz^ .i,t^\]zi 



+h.c. 



9 



8(27r)' 



+ (+ 



rTr 



(5.3) 



+h.c. 



+ (+ 



,<-} + ^fe,v;r}) + [[z'+,z'+],z+]) 

■ -) (5.4) 
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Note that the last term of the second hne in (|5.3| ) is just the part for the subsector 
su(2|3) discussed in [17]. 

We only consider planar diagrams in this paper, so the terms we are interested in must 
have adjacent fermions and bosons, and only has non-zero correlation function with 

Of~^ or Of~^ at planar level to 0{g^). The matrix model vertices must involve all the 
impurities in both operators in the two-point correlation function. So, using the cyclic prop- 
erty of the trace of operator string, for effective vertices like Tv{'ij;^ ij;^(j)f (p'j' Z~^) in which 
the two impurity bosonic fields are adjacent, we can change the position of Z or Z from the 
end of the three bosonic fields to the beginning of them or vice versa without changing the 
two point correlation function, that is il)^(j)f (f)'^ Z^) Tv{il)'j^ilj^^Z^(j)f(f)'^). Note 

that this operation only works for two impurity operators. By appropriately using this 
operation, we collect the terms we are interested in as: 

HiA) + H^B) - TTf^Trj - ( ^ <t>i + 2Z- Z' zA - ^[^"+^'^1+) 

' I ^«=l, 2,4,5 ^ 

J,a+ 



-I 



-I 



(2<A-05-Z- + 2cp-^cp-^Z- - cp-^Z-cp- - cp-^Z~cp-y{ptJt+ - i^+Ji+) 

+(/i.c.,+ ^ -)| , (5.5) 

where 4'[i4'j] = \{'Pi4>j ~ 'Pj'Pi)- We can define: 

/nP'-' — roP'-^ _ dP'-^ dP'-^ — roP'-^ -l nP''^ 

= If^'if + ^^li) ' ^Ei] = li^'lf - K') ■ (5-6) 

From ( ^I^ ) we see that O^'"^ are mixed with ^'j'^^'g] ^ and ^'j'^^j ^, Ops' 14] ™ixed with ^'^{3) ^' 

^(12 45) with ^('{;^22)' ^® ^^^^ these mixed operators in the rows of table |. 

These mixing patterns are consistent with group analysis as follows: if we keep Z = 
2$34 and Z = 2$i2 invariant, the SU{A)r symmetry is broken to 5'0(4) = SU{2) x SU(2) 
with '0^'^ and -(/ja 4 as fundamental representations of the two S'C/(2)'s respectively. Using 



(A. 8) and omitting Z, we can obtain: 

^C'[14,25] ~ <^^'^^^^ - «>23$14 ^ $24^13 _ ^13^24 ^ ^5 7^ 

-0(12,45) ~ ^'''J''' - ^'''^''^ - + , (5.8) 

-zO(24,i5) ~ 'J'''^''' + '^^^^^^ - ^^""^^^ - ^^^^^^ , (5.9) 
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^[121 ^[341 


*'^[14,25l 


^(12) 


(12,45) 


^— — T 

(11,22) 


^'-^(24,15) 


^P,J-1 
[11,22] 


''-^(14,25) 


^(34) 


fl2,45l 


^(33,44) 


''-^[15,241 


^P,J-i 
[33,44] 



Table 1: The operators on the same row will be mixed. 



from which we see that the indices 3 and 4 are antisymmetric and form a singlet under 
the second SU{2). The indices 1 and 2 have precisely the same transformation property 
as the corresponding mixed operators in table |l|. The same analysis can be performed on 
the other rows. In fact, ^'n, ^'22 and ^'(12) form a 3 representation of 50(4) and ^'33, ^'44 
and ^'(34) form a 3 representation. The corresponding scalars are also expected to form 3 
and 3 representions. As a result, the operators V^'"^ are only mixed with the first line of 
table but not with the other lines. 

6. Mixing of operators 



Let us look at the mixing between Oi, ^[12] and ^'[34] first. Since there is no correlation 



[34] 



between \I'[i2] and ^[34] up to 0{g^), we can consider the mixing between Oi and ^'[12], and 
Oi and ^'[34] separately. Since the coefficients before ^^[1^2] and ^/^l^^^l in the vertices differ 
only by a minus sign, the mixing coefficients for ^[12] and ^'[34] with Oi must differ only by 
a minus sign. For operator basis {Oi^''\Oi''\ . . . ,0^^'^^''^ ^ ■ ■ ■ ^"^[12]^ ^'^ the 

H matrix should combine the one-loop H matrix for O^'"^ and ^|']^2]^^ and some 0{g^) 
crossing terms: 



for J odd: 



9^ 
47r2 



-1 2 -1 
-1 1 



x/2(27r) ^2(2 



V 



2y2(27r) 

g^/lv 

2v^(27r) 



\ 



-1 2 -1 
-1 2 



(6.1) 



After diagonalization, we require the non-diagonal elements to be 0{g^). So we require 
the correction to the eigenvector Oi to be 0{g). The upper right non-diagonal element 
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contribution to the eigenvalue is aheady 0{g'^). The upper left block of the H matrix 
acting on the 0(1) part gives just the one-loop case which determines the eigenvalue to be 
(g'^A^/47r^)4sin^(n7r/(J + 3)). As stated at the end of section ^, we could make the ansatz 
for the coefficients before ^^l2]~^ *° 2gbsm !2P+^hlL_ xhe boundary condition from the 
last row of the H matrix then determines a = 4. The eigenvalue equation for the coefficient 
before ^!^'-^~^ acts as another boundary condition which determines b. We then obtain the 
eigenvector with 0{g) correction to O^'"^: 

^J,n .?,J,n , V2gVN . nvr . (2p + 4)ri7r j_i 

p=0 

for n = 0, . . . , [J/2] + 1 . (6.2) 



We can also require the coefficients before ^'|']^2]~^ *° and those before O/ to be 

0{g). This time the lower left nondiagonal element contribution is already 0{g'^), and the 
upper right nondiagonal elements contribute to two eigenvalue equations. As a result, the 
coefficient before O^^' should be considered as independent of the other coefficients: we 
set it to be A^i and for p = 0, . . . , (J — l)/2 the coefficients before O^'^ are 26cos (?P+^l!lZL, 
As in previous case, the last line in the upper left block determines o = 1. Those two 
equations from the first two lines determine A-i to be zero and the value of b, thus the 
eigenvector in this case is 

^J-i,n ,T,J-i,n gVN . nvr {2p+l)mT pj 
^r-|„, ' = Wr-,„, ' — — -sm— -> COS ; Cr , 

forn = 1, . . . , [( J + l)/2] . (6.3) 



It is easy to check that for the J = even case the eigenvectors can also be written in this 
form. Notice that the total number of eigenvectors in ( |6.2| ) and ( |6.3D is equal to the rank 
of the matrix and the eigenvalues are not degenerate. So these are all the eigenvectors of 
the H matrix (|6.1| ). 

For the mixing with ^'[34], only the nondiagonal blocks in ( |6.1| ) change signs. As a 
result the coefficients of 0{g) of the eigenvector change signs. So combined with (|6^), we 
obtain the 0{g) correction to the eigenvector O^' 

^J,n .?.J,n , V2g\fN . mr . (2p + 4)ri7r / j_i ,.,p,j-i\ .^ 

Qi =0, +^^^-7^L^-— 7^^(^[i2] -^[34] j (6-4) 

p=0 

This is consistent with equation (3.17) in [20]. ( |6.3D is consistent with the first two terms 
in equation (2.12) in [26]. To obtain the full 0{g) correction, one must also consider the 
mixing matrix elements between ^'[12] and D Z"^ Z Z'^ . 

As for the mixing between Of''^„^. and ^^{i)^^ have the H matrix for 
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operators basis (^^£'['1425]' 
for J odd: 

/ 2 



,iO 



[(J-1)/2],J ^0,J-1 



[14,25] 



(12) 



9^ 

47r2 



^[(J-l)/2],J-l 
' (12) 



2v^{27r) 



(6.5) 



-1 2 -1 
-1 3 



V2gVN __gVN_ 
(27r) V2(27r) 



V 



-1 2 -2 
-1 2 



(6.6) 



It is easy to obtain the eigenvectors as before: 
J 

'[14,25] ~ " 



^ J,n 
^(12) 



p=0 

J,n 



J + 2 



[^^'25] ^(2^) 
^TT „o,J+l 



g\/ N nvr 
sin 



i4 (2p + 3)n7r j_i 

jT2>^^°^ j + 2 ^^a^) 

p=0 



nvr 



J 



H ^ ( — cos — Om^i'Tki + sin — y sin ^ ^ J — z — Ofi 

(^2) V2(2^)V "^ + 1 ' ' ^' '^ + ^^ 



(2p + 2)n7r j+i 



J + 2 



[14,25] 



(6.7) 



As in previous case the coefficient before the O^''^^^ term should be considered indepen- 
dently. This term corresponds to the non-asymptotic term in (3.1) of [19]. The other 
operator mixing cases in table || can be easily obtained by replacing ^^'^ ^('12) 

by the corresponding operators in the same column in table 0, because they produce the 
same H matrix. As we mentioned in the last section, these operators are not mixed with 
Df,ZZ ...D^'ZZ ...Z. 



7. Discussion 

We have discussed the mixing between single trace BMN operators in the pure scalar sector 
and operators with one scalar fermion pair. A kind of matrix model effective vertices up to 
0{g^) are deduced during the discussion. These effective vertices can also be used to discuss 
non-planar diagrams. All the mixing patterns at planar level are reflected in the vertices 
( |5.5| ). The 0{g^) matrix elements in the eigenvalue problem act as a boundary condition 
to the 0{g'^) matrix. The results are consistent with [20] and [26]. It is worth pointing out 
that even though the correlation functions between operators with nonadjacent fermions 
and pure scalar operators are zero up to 0{g^) at planar level, they are mixed because the 
mixings between operators with adjacent fermions and those with nonadjacent fermions 
are needed to diagonalize the one-loop dilatation operator. 
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Since we discuss only operators with one fermion pair, the 0{g^) effective vertices do 
not include the interactions between two fermion pairs like: 



A A 
/ \ I \ 

\ H \ 



If one wants to discuss operators with two or more scalar fermion pairs, one has to include 
the terms from this diagram which can be calculated as in [9]. In general, this term could 
also be obtained by using the results in [27]. The 0{g^) terms ( |5.3| ) and ( |5.4| ) can be 
directly used in more than two fermion pair cases. It is an interesting problem to discuss 
the integrability of this kind of dynamic spin chain to 0{g^). Unlike the 5'C/(2|3) spin 
chain in [17,19,29], the fermions are not invariant under the SU{'i) subgroup of 5'C/(4)/j. 
To be precise, they transform under SU{2) x SU{2) subgroup. Understanding how to use 
the coordinate space Bethe ansatz in this kind of dynamic spin chain could be a further 
direction for research. 

One can also discuss the mixing between the SO{A) singlet parts of '^^'^ and V'^''^ 
along these lines. This requires calculating the two point correlation function to 0{g^) and 
we could expect that the 0{g^) matrix elements act as boundary conditions, just as in the 
cases we discussed here. One can also go further, calculating the dilatation operator matrix 
between O^'^ and V^'^ to 0{g^), and we still expect that the 0{g^) matrix elements act 
as boundary conditions. These will produce the full equations (2.1) and (2.12) in [26]. 

A few remarks on the method in [20]: in [20], the authors use the supercurrents 
deduced from dimensional reduction of the D = 10 SUSY YM. But if one discusses the 
quantum properties of = 4 SUSY YM, one must include the gauge fixing and ghosts in 
the quantized Lagrangian. One would expect that this should modify the supercurrents, 
which could be different from those in [20] in the gauge sector and ghost sector. The 
problem is more explicit in = 1 SUSY YM theory. We can do gauge fixing and add 
ghosts in superfield formalism without breaking M = 1 supersymmetry. Then we can 
have a supercurrent involving gauge fixing and ghosts. On the other hand, for the original 
lagrangian without ghosts and gauge fixing, we also have a supercurrent. Which one should 
be chosen to calculate the higher order correction to the SUSY variation of a composite 
operator? The method in [20] is equivalent to choosing the supercurrent without ghosts. 
Since all the perturbation calculations are based on the Lagrangian after quantization, 
the more consistent choice is to use the former supercurrent. It turns out that in [20] 
when they only discuss the scalar and fermionic sectors, using those supercurrents from 
ten-dimensional SUSY YM is fine. We can see intuitively why this is not a problem. After 
careful scrutiny, we can conclude that their method is consistent using either of these two 
currents. This is because after the Faddeev-Popov quantization the partition function is 
essentially proportional to the path integral using the classical action divided by the volume 
integration of the gauge redundancy: 

m.,) . . . OK-0>.-P = j[n.,^.A.(,)],[A(,)] • (^-1) 
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This is equation (15.5.20) of [30]. Operators Oi are supposed to be gauge invariant. The 
subscript F-P means the path integral after Faddeev-Popov quantization and subscript c 
means the classical one. p[A(x)] is the invariant (Haar) measure on the space of group 
parameters, (pn here denotes generally all the fields to be integrated. Note that the de- 
nominator and the constant C are gauge invariant. We could deduce the Ward Identity 
using this path integral with classical action and get back to Faddeev-Popov ghosts later. 
Then the current in the Ward Identity is the classical one without ghosts: 



d^{j>!{x)Oi{xi)...Oiixi))^_^ 

_ C/[n„,,#„(x)]5^j,^(x)Oi(xi) . . . Oi{xi)e'^'^ 



f[U^,,A^ix)]p[A] 
C /[n„,,.#„(x)]Oi(xi) . . . Oi{xi)5r{x)j^^ 



5S ^iSc 



/[n„,,A"(x)]p[A] 
C J[Un,,dMx)]Oi{xi) . . . SO,{x,)6Hx - X,) . . . Oiixi)e'S^ 



/[n„,.A"(x)]p[A] 
; Y,{Oi{xi) . . . SOi{xi)6\x -Xi)... Oi{xi)), 



(7.2) 



Notice that the infinite integral in the denominator and numerator in (7^) should be 
understood as being put on a finite lattice and then taking the infinite limit in front of the 
fraction. Since the Faddeev-Popov integral is well-defined, the limit is well-defined. This 
can be done in every line of ( [7. 2] ). In their method, the Ward Identity is used on some 
operators which are not gauge invariant to deduce the SUSY transformation property of 
one operator. Then this SUSY transformation is used to obtain the SUSY variations of 
some gauge invariant operators. This is equivalent to using the Ward Identity on the gauge 
invariant operators from the beginning, thus (|7.2|) can be applied. 
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A. TV = 4 SYM convention 

The conventions we use in this paper follow closely the paper [20] except for some normal- 
ization differences: 

Tr(T'^r^) = 5''^ , [T^^*] = V2if^^T^ , (r'');(T'^)f = 6\6^j (A.l) 
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We choose Minkowski metric Qui, = diag{+, —, —, —). The M = A Yang-Mills La- 
grangian is 



L = Tr 



where 



There are some useful relations: 



The propagators are: 



(A.2) 



2' 

For scalar s, we use Zi, and in different contexts: 

^ _? U -'7 for z,i,fe = 1,2,3 (A.3) 

— 2^jkiZ I 5>i4 — 2 A 



Zi = ^i(j)i+ i(f)i+3) , = -^(</)i - i(f>i+3) . (A.4) 



+ Z^Z'^ = <l)t<l>) + C+30'+3 (A.5) 

Zazbj _ z-iz] = -i{<l>t<l>]+z - (A.6) 

4 + Z-'z'i = <t>1<t>] - </.f+30'+3 (A.7) 

4 - Z-Z"^- = i{ct>U]+^ + ^r+a*/-?) (A.8) 



{Z'^{x)Z\y)) =5'^''l^^y, 

l,Z\x)Z\y)) = l,Z\x)Z\y)) = , (A.9) 



where 
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